NAVAL  POSTGRADUATE  SCHOOL 
Monterey,  California 


THESIS 


DESIGN  OF  A  TOROIDAL  THERMOACOUSTIC 
PRIME  MOVER 

by 

Yang,  Szu-Wei 

June  1995 

Thesis  Advisor; 

Thesis  Co-Advisor: 

Anthony  A.Atchley 
Thomas  J.  Hofler 

Approved  for  public  release;  distribution  is  unlimited 

.Cl  ^ 


19960122  105 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  0704 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for 
reviewing  instruction,  searching  existing  data  sources,  gathering  and  maintaining  the  data  needed,  and  completing  and 
reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this 
collection  of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  headquarters  Services,  Directorate 
for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  222024302,  and  to  the 
Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188)  Washington  DC  20503. _ 

I.  AGENCY  USE  ONLY  flMiw  */anifc>  2.  REPORT  DATE  3.  REPORT  TYPE  AND  DATES  COVERED 

June  1995  Master’s  Thesis _ 

4.  TITLE  AND  SUBTITLE  5.  FUNDING  NUMBERS 

DESIGN  OF  A  TOROIDAL  THERMOACOUSTIC  PRIME  MOVER _ 

6.  AUTHOR(S) 

Yang,  Szu-Wei _ 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES)  8.  PERFORMING  ORGANIZATION 

Naval  Postgraduate  School 

Monterey  CA  93943-5000 _ REPORT  NUMBER _ 

9.  SPONSORING/MONTTORING  AGENCY  NAME(S)  AND  ADDR£SS(ES)  10.  SPONSORING/MONITORING 

AGENCY  REPORT  NUMBER 

II.  SUPPLEMENTARY  NOTES 

The  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect  the  official  policy  or  position  of  the 
Department  of  Defense  or  the  U.S.  Government. _ _ 

12a.  DISTRIBUTION/AVAILABILITY  STATEMENT  12b.  DISTRIBUTION  CODE 

Approved  for  public  release; 

Distribution  unlimited. _ 

13.  ABSTRACT  (maximum  200  words) 

The  topic  of  this  thesis  is  the  design  of  a  toroidal  prime  mover.  The  primary  design  tool  is 
a  MATLAB  program  based  on  a  standing  wave  analysis  of  thermoacoustics.  A  conventional  prime 
mover  that  has  been  well  studied,  was  used  for  validating  the  program.  DeltaE,  a  program  developed 
at  Los  Alamos  National  Laboratory  for  designing  thermoacoustic  engines,  was  also  used  as  an 
auxiliary  tool  for  validation  test  purposes.  The  validation  suggested  that  the  MATLAB  program 
models  the  onset  conditions  of  a  conventional  prime  mover  with  fair  accuracy.  The  program  was  then 
applied  to  a  toroidal  prime  mover.  A  series  of  analyses  were  conducted  on  a  toroidal  prime  mover  to 
determine  the  influence  of  prime  mover  circumference,  stack  length,  and  duct  radius.  The  final  design 
calls  for  a  150  cm  circumference,  2.54  cm  duct  radius,  and  1.9  cm  stack  length.  The  contributions  of 
this  thesis  are:  (1)  to  develop  and  validate  the  MATLAB  simulation  program;  (2)  to  provide  a  draft 
design  for  a  toroidal  prime  mover  for  future  research. _ 

14.  SUBJECT  TERMS  15.  NUMBER  OF  PAGES 

Thermoacoustics;  Thermoacoustic  Heat  Transport;  Prime  Mover;  Toroidal  Prime  Mover  62 

16.  PRICE  CODE 

17.  SECURITY  18.  SECURITY  19.  SECURITY  20.  LIMITATION  OF 

CLASSIFICATION  OF  CLASSIFICATION  OF  CLASSIFICATION  OF  ABSTRACT 

REPORT  Tins  PAGE  ABSTRACT  UL 

Unclassified  Unclassified  Unclassified _ _ 


NSN  7540-01-280-5500 


Standard  Form  298  (Rev.  2-89) 
Prescribed  by  ANSI  Std.  239-18 


Approved  for  public  release;  distribution  is  unlimited. 

DESIGN  OF  A  TOROIDAL  THERMO  ACOUSTIC 
PRIME  MOVER 


Yang,  Szu-Wei 

Lieutenant  Commander,  Taiwan  Navy 
B.S.,  Chinese  Naval  Academy,  1983 

Submitted  in  partial  fulfillment 
of  the  requirements  for  the  degree  of 

MASTER  OF  SCIENCE  IN  ENGINEERING  ACOUSTICS 

from  the 


NAVAL  POSTGRADUATE  SCHOOL 
June  1995 


iii 


IV 


ABSTRACT 


The  topic  of  this  thesis  is  the  design  of  a  toroidal  prime  mover.  The  primary 
design  tool  is  a  MATLAB  program  based  on  a  standing  wave  analysis  of  thermoacoustics. 
A  conventional  prime  mover  that  has  been  well  studied,  was  used  for  validating  the 
program.  DeltaE,  a  program  developed  at  Los  Alamos  National  Laboratory  for  designing 
thermoacoustic  engines,  was  also  used  as  an  auxiliary  tool  for  validation  test  purposes. 
The  validation  suggested  that  the  MATLAB  program  models  the  onset  conditions  of  a 
conventional  prime  mover  with  fair  accuracy.  The  program  was  then  applied  to  a  toroidal 
prime  mover.  A  series  of  analyses  were  conducted  on  a  toroidal  prime  mover  to 
determine  the  influence  of  prime  mover  circumference,  stack  length,  and  duct  radius.  The 
final  design  calls  for  a  150  cm  circumference,  2.54  cm  duct  radius,  and  1.9  cm  stack 
length.  The  contributions  of  this  thesis  are:  (1)  to  develop  and  validate  the  MATLAB 
simulation  program;  (2)  to  provide  a  draft  design  for  a  toroidal  prime  mover  for  future 
research. 
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I.  INTRODUCTION 


There  have  been  many  NFS  theses  involving  several  aspects  of  thermoacoustic  heat 
transport.  One  area  of  particular  concentration  has  been  thermoacoustic  prime  movers 
[See  Refs.  1-3,  for  example].  A  thermoacoustic  prime  mover  is  a  class  of  thermoacoustic 
heat  engine  that  converts  thermal  energy  into  useful  work  in  the  form  of  sound.  The 
purpose  of  this  thesis  is  to  explore  the  design  of  a  toroidal  prime  mover. 

A  typical  prime  mover  configuration  is  illustrated  in  Fig.  1.  During  operation,  a 
temperature  difference  is  imposed  between  the  two  heat  exchangers,  resulting  in  a 
temperature  gradient  in  the  stack.  (The  terms  hot  and  cold  are  used  to  indicate  that  there 
is  a  temperature  difference  between  these  regions  of  the  prime  mover.  They  are  not  meant 
to  imply  anything  about  the  absolute  temperature  involved).  When  the  gradient  is 
sufficiently  high,  acoustic  oscillations  spontaneously  build  up  in  the  prime  mover.  At  this 
point  the  prime  mover  is  said  to  have  reached  onset.  The  conditions  necessary  for  a  prime 
mover  to  reach  ortset  are  typically  discussed  in  terms  of  which  of  two  quantities  is  greater: 
the  thermoacoustic  gain  of  the  stack,  or  the  dissipation  due  to  thermo-viscous  wall  losses 
and  power  delivery  to  an  external  load.  The  importance  of  each  of  these  factors  is 
determined  by  the  geometry  of  the  prime  mover,  e.g.,  where  the  stack  is  located  relative  to 
the  nearest  acoustic  pressure  antinode,  etc.  The  positions  of  the  components  of  a  prime 
mover  relative  to  the  acoustic  field  are  usually  fixed  by  incorporating  well  defined  acoustic 
boundary  conditions  into  the  design  of  the  engine.  For  instance,  placement  of  a  rigid 
boundary  in  the  engine  insures  that  the  acoustic  velocity  is  zero  at  that  point.  These 
"built-in"  boundary  conditions  serve  as  convenient  starting  points  for  computational 
analysis  of  the  performance  of  an  engine. 

A  frequently  asked  question,  when  explaining  how  a  prime  mover  generates  sound 
(once  the  temperature  gradient  in  the  stack  is  sufficiently  high)  is  "What  is  the  origin  of  the 
low  amplitude  sound  that  gets  amplified  by  the  prime  mover?"  When  answering,  one 
usually  resorts  to  postulation  the  existence  of  fluctuations  or  noise  having  a  component 
with  the  right  frequency  and  acoustic  phase  to  be  amplified  by  the  stack. 
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Figure  1  -  A  typical  prime  mover  configuration 
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In  light  of  this  discussion,  a  toroidal  prime  mover  poses  some  interesting  questions. 
As  shown  in  Fig.  2,  a  toroidal  prime  mover  has  the  cold  end  connected  to  the  hot  end  via  a 
duct  that  does  not  impose,  at  first  glance,  any  obvious  dominating  boundary  conditions 
that  would  force  an  antinode  at  any  particular  point.  Indeed,  if  the  stack  and  heat 
exchangers  were  not  there  and  if  the  toroid  had  a  uniform  cross  section,  a  standing  wave 
should  have  no  preferential  orientation.  In  other  words,  there  is  no  place  that  one  knows 
a  priori  is  going  to  be  a  pressure  antinode,  for  example.  However,  one  might  argue  that 
even  with  periodic  boundary  conditions,  fluctuations  still  exist  and  there  will  be  a 
component  that  matches  the  conditions  necessary  for  amplification. 

Assuming  that  a  toroidal  prime  mover  can  go  into  onset,  a  number  of  interesting 
questions  arise.  Will  the  standing  wave  (if  it  is  an  standing  wave)  adjust  its  spatial 
orientation  to  optimize  the  thermoacoustic  gain?  In  other  words,  will  the  prime  mover  put 
the  stack  where  it  wants,  where  its  performance  is  optimized?  A  rigid  ended  prime  mover 
does  not  offer  this  possibility  of  self-optimization.  Will  deviations  in  measured  stack 
locations  from  those  predicted  be  indicative  of  unaccounted  for  processes  such  as  acoustic 
streaming  or  turbulence?  Another  interesting  question  concerns  multimode  excitation.  It 
has  been  previously  shown  that  two  modes  can  be  excited  in  a  prime  mover 
simultaneously  [Refs.  3  and  4].  However,  due  to  the  presence  of  rigid  boundaries,  the 
stack  placement  is  optimized  for  only  one  mode.  A  toroidal  geometry  with  periodic 
boundary  conditions  offers  the  possibility  that  each  mode  could  optimize  itself 

All  these  questions  presuppose  that  the  prime  mover  will  go  into  onset.  The 
presence  of  the  stack  will  break  the  ambiguity  of  the  standing  wave's  spatial  orientation. 
The  stack  has  the  highest  flow  impedance  in  the  toroid.  Because  this  impedance  wilt  tend 
to  favor  low  acoustic  velocities,  the  perturbation  due  to  the  presence  of  the  stack  might 
force  a  velocity  node  at  the  stack.  This  condition  would  necessarily  make  any 
thermoacoustic  gain  zero,  preventing  the  prime  mover  from  ever  going  into  onset. 
Further,  the  fundamental  mode  of  a  uniform  cross  section  toroid  (one  without  a  stack), 
corresponds  to  the  circumference  being  equal  to  one  acoustic  wavelength.  There  are  two 
orthogonal  degenerate  modes  that  satisfy  this  condition.  The  presence  of  a  nonuniformity 
will  break  this  degeneracy.  One  mode  will  have  a  velocity  node  at  the  stack,  the  other  will 
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have  a  pressure  node  at  the  stack.  Neither  of  these  modes  will  be  supported  by 
thermoacoustics,  but  linear  combinations  would  be  supported.  Other  nonuniformities  in 
cross  section,  will  also  have  an  effect  on  the  spatial  distribution  of  the  acoustic  variables. 
Hence,  by  placing  a  constriction  or  expansion  in  the  toroid,  the  previously  mentioned 
velocity  and  pressure  nodes  will  be  displaced  from  the  stack,  allowing  for  the  possibility  of 
thermoacoustic  gain. 

This  discussion  leads  one  to  conclude  that  designing  a  functioning  toroidal  prime 
mover  may  be  considerably  more  difficult  than  designing  a  conventional  prime  mover. 
The  research  described  in  this  thesis  addresses  only  the  first  part  of  the  problem.  That 
problem  is  to  design  a  toroidal  prime  mover  that,  presupposing  onset  is  possible,  should 
reach  onset  at  a  temperature  difference  that  is  easily  attainable  in  the  laboratory. 

The  approach  taken  is  to  use  a  simple  formulation  of  the  basic  equations  for 
thermoacoustics  to  predict  the  onset  conditions  for  a  toroidal  prime  mover.  The  program 
developed  around  this  formulation  \vill  first  be  validated  by  comparing  its  predictions  to 
data  taken  previously  on  a  rigid  ended  prime  mover  [Ref  1].  The  predicted  onset 
conditions  will  also  be  compared  to  those  determined  using  DeltaE,  a  program  developed 
at  Los  Alamos  National  Laboratory  for  designing  thermoacoustic  engines.  After 
validation,  the  MATLAB  program  will  be  used  to  design  a  toroidal  prime  mover. 

Chapter  II  contains  a  brief  review  of  the  theory  pertinent  for  this  work.  The 
validation  of  the  program  is  discussed  in  Chapter  III.  Some  experiments  with  a  toroidal 
prime  mover  are  presented  in  Chapter  IV.  An  improved  design  is  presented  in  Chapter  V. 
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n.  THEORY 


It  has  been  shown  previously  that  tracking  the  quality  factor  (Q)  of  a  prime  mover 
as  a  function  of  temperature  difference  across  the  stack  is  a  convenient  way  to  predict 
onset.  This  chapter  contains  a  discussion  of  how  to  calculate  the  Q  of  a  prime  mover  as  a 
function  of  the  applied  temperature  difference.  The  discussion  is  based  on  a  standing 
wave  analysis  described  in  detail  in  Ref  5  and  draws  heavily  upon  the  derivation  presented 
there. 

In  addition,  brief  descriptions  of  the  two  different  computer  programs,  one  written 
in  MATLAB  [Ref  6]  by  the  author,  and  DeltaE  [Ref  7]  will  be  given  at  the  end  of  this 
chapter. 

A.  STANDING  WAVE  ANALYSIS  OF  A  THERMOACOUSTIC  PRIME 
MOVER 

The  derivation  begins  by  making  the  connection  between  the  power  output  of  a 
prime  mover  and  quality  factor.  Next,  an  expression  is  derived  for  the  power  output  of  a 
prime  mover  stack.  Following  this  derivation,  the  thermal  and  viscous  losses  in  the 
remainder  of  the  prime  mover  are  discussed.  Finally,  an  expression  is  derived  for  the 
energy  stored  in  the  prime  mover. 

1.  Relationship  Between  Power  Output  And  Quality  Factor 

The  reciprocal  of  the  quality  factor  of  a  prime  mover  below  onset  can  be  expressed 
as 


\IQ  =  -WI(iiEsT 


(1) 


The  net  power  output  W  is  the  sum  of  the  power  output  of  the  individual  sections 
of  the  prime  mover : 


W=Who,  +Wt,othx  +Wp„,s  +Wcoldhx  ^Wcold  .  (2) 

The  subscripts  hot,  hot  hx,  pms,  cold  hx,  and  cold  refer  to  the  hot  end,  the  hot  heat 
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exchanger,  the  prime  mover  stack,  the  cold  heat  exchanger  and  the  cold  end,  respectively. 

Est  is  the  acoustic  energy  stored  in  the  prime  mover,  w  is  the  angular  frequency  of  the 
acoustic  oscillation.  At  onset  fV=  0,  so  l/Q  =  0.  W=  0  means  that  the  power  output  of 

the  stack  exactly  cancels  the  losses  elsewhere  in  the  prime  mover,  so  the  net  power  output 
is  zero. 

2 .  Acoustic  Power  Output  Of  The  Prime  Mover  Stack 

In  Eq.  (1)  fVp,„s  is  the  power  output  of  the  stack.  In  general,  fFp,„s  fVp,„s , 

where  d  fVp,„s  is  the  differential  power  output  of  a  segment  of  the  stack  located  between  x 
and  x  +  dx.  It  is  found  by  taking  the  difference  in  the  acoustic  intensity  at  the  two  ends  of 
the  segment  and  multiplying  the  result  by  the  cross-sectional  area  of  the  fluid  in  the  stack 
‘S’pmi(the  X  direction  is  along  the  resonator  axis,  the  direction  of  sound  propagation). 

^  ^ oms  ~  S pms\^U^ x+dx  ~  (3) 

The  overbar  indicates  time  average  and  the  brackets  <)  indicate  average  in  the  y  direction. 
(The  y  direction  is  perpendicular  to  the  resonator  axis  and  the  plate  surface.)  The  stack 

plates  are  2/  thick  and  spaced  by  2yo.  Expanding  (pu}x+dx  and  assuming  that  p  is 
independent  of^',  Eq.  (3)  can  be  written  as 


dWp^s  =Sp„s^^dx.  (4) 

The  acoustic  variables  can  be  expressed  as p  =pie^  and  u  =  uie^,  where 

pi  =P4C0s[(t)(x)],  (5) 

and 


In  these  equations.  Pa  is  defined  to  be  the  peak  acoustic  pressure  amplitude  at  the  rigid 
end,  at  jc  =  0.  Changes  in  pressure  amplitude  due  to  changes  in  resonator  cross  section  are 
neglected.  Also,  (|)(x:)  represents  the  phase  of  the  standing  wave  at  the  location  of  the 
segment.  If  the  propagation  constant  k  is  independent  of  position,  (t>(x)  =  kx.  However,  in 
general. 


<|)  W  =  lo 


(7) 
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The  term  (1  +  Hyo)  arises  from  the  requirement  that  volume  velocity  be  conserved  and  the 
desire  to  relate  the  acoustic  velocity  in  the  stack  to  Pa  ■ 

The  time  average  of  the  product  of  two  complex  quantities  such  as  pand  <w>  is 


p{n)  =  \Re\px{ui)] 


(  ~  denotes  complex  conjugate),  Eq.  (4)  becomes 


dWp„,,  =\Sp^sRe{px^^^^{ux))dx 


The  expression  for  dpxidx  and  d(px)ldx,  are  obtained  from  Swift  [Ref  8]: 


dpx  -mm,  , 

=  “i — 
dx  1  -/v 


d{ux)  _  j  fk-fv  VT„dpx  j  0)2  r  Y-l.V 

(Jx  (1 -a)(l +8j)  Tm  dx  1+e/v^*' 


The  functions  fk,fv  and  Sj  are  given  by 


tanh[(l  +y)yo/5t] 

(1  +y)yo/5ifc 

tanh[(l  +y);^o/5v] 

(1  +y>o/6v 

p^Cp6ifctanh[(l  +7)>^o/5t] 

p,C55itanh[(l +y)//5s] 


where  the  thermal  and  viscous  penetration  depths  in  the  gas  are  defined  as 


5*=  J2K/p„Cp(}i , 
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6v  =  y2|Li/p,„© . 


(16) 


The  definition  of  the  thermal  penetration  depth  in  the  solid  plate  5s  is  similar  to  that  of  6*. 
The  Prandtl  number  c  is  given  by 


a  =  Cpp/AT, 


(17) 


Substituting  Eqs.  (5),  (6),  (10),  and  (11)  into  Eq.  (9)  gives  the  power  output  of  the 


element 


^  ^pms  —  ^SpmsR^ 


I  Pjo) 


[  PmCl^ 
a 


Jk-fv  \+llyo  VT,n^  ...  ,  (Y-1)~ 

(l-o)(l+g,)A:r„  '  1+8,-^^ 


xcos2(x)  -yil~^sin^(l)(jc) 

^  /v 


^dx. 


(18) 


The  total  acoustic  power  output  of  the  stack  is  found  by  integrating  Eq.  (18)  along  the 
entire  length  of  the  stack. 

3.  Thermal  and  Viscous  Losses  In  The  Remainder  of  The 
Prime  Mover 

Thermal  and  viscous  losses  can  be  calculated  from  Eq.  (18)  with  only  slight 

modifications.  The  calculation  is  simplified  by  the  absence  of  temperature  gradients  in  the 
remainder  of  the  prime  mover,  i.e.  Vr,„  =  0. 


=\Sh.-cRe\^^ 

2  I  p,„a^ 


^  j  sin^(t)(r) 


^dx.  (19) 


Equation  (19)  also  describes  the  losses  in  the  cold  end  and  hot  end  if  /  is  set  to  zero  and 
the  appropriate  form  for  S,fk,  and/v  are  used.  For  a  circular  tube  of  radius  R, 

2Jx[(J-\)RI5k\ 

(/-I  )(/?/6*>/o  [(/■  -  1  )/?/5t]  •  (20) 

and 
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(/-IKTJ/SvVolO-lWSv]' 


(21) 


The  thermal  losses  at  the  end  caps  are  found  by  evaluating  Eq.  (19)  at  x=0  (the  cold  end  ) 
and  x=L  (the  hot  end)  and  replacing  >^0  with  L  . 

4 .  Stored  energy 

The  energy  stored  in  the  prime  mover  is  found  by  integrating  the  time-averaged 
acoustic  energy  density  throughout  the  entire  volume  of  the  prime  mover.  The  result  is 


where  Xcoid  and  Xhot  refer  to  the  length  of  the  cold  and  hot  end,  respectively. 

B.  MATLAB  PROGRAM 

A  MATLAB  program  was  written  to  implement  the  standing  wave  analysis 
presented  above  to  calculate  \/Q  vs.  AT.  In  regions  where  there  is  no  gradient 
(everywhere  except  in  the  stack),  integrals  are  evaluated  analytically.  A  linear  gradient  is 

assumed  in  the  stack.  The  stack  is  divided  into  10  isothermal  segments.  The 
temperature-dependent  thermal-physical  parameters  and  acoustic  phase  are  evaluated  for 
each  segment  and  the  work  output  computed.  The  total  work  output  of  the  stack  is  the 
sum  of  the  work  output  of  each  segment.  In  rigid-ended  prime  movers,  the  position  of  the 
elements  of  the  engine  relative  to  the  standing  wave  are  determined  by  boundary 
conditions.  In  the  toroidal  prime  mover,  the  location  of  the  elements  relative  to  the 
standing  wave  is  an  adjustable  parameter. 
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C.  DELTAE 


DeltaE  (Design  Environment  for  Low  Amplitude  Thermoacoustic  Engines)  is  a 
computer  program  written  by  Ward  and  Swift  at  Los  Alamos  National  Laboratory  [Ref  7], 
It  solves  Rott's  wave  equation  for  a  geometry  defined  by  the  user.  More  information  on 
DeltaE  can  be  found  on  internet  at  http://rott.esa.lanl.gov. 
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m.  VALIDATION 


The  goal  of  this  thesis  is  to  design  a  toroidal  thermoacoustic  prime  mover.  As 
discussed  in  the  previous  chapter,  the  primary  design  tool  is  a  MATLAB  program  based 
on  the  analysis  presented  in  Ref  5.  The  primary  design  goals  are  1)  to  have  an  onset 
temperature  difference  that  is  easily  reached  experimentally  and  2)  to  have  dimensions  that 
are  practical  from  machining  and  experimental  considerations. 

The  purpose  of  this  chapter  is  to  explain  the  process  used  to  validate  the  MATLAB 
program.  The  program  was  validated  in  two  ways.  1)  It  was  used  to  predict  the 

performance  of  Lin's  prime  mover  below  onset,  a  system  that  has  been  well  studied 
previously.  2)  Predictions  of  onset  temperature  difference  (ATonset)  were  compared  to  the 
predictions  of  DeltaE. 

A.  LIN’S  PRIME  MOVER 

This  prime  mover,  described  in  detail  in  Ref  1  and  Ref  5,  was  used  originally  to 
investigate  the  low  amplitude  work  output  of  a  prime  mover  stack  as  a  function  of  mean 
gas  pressure  and  temperature  difference  across  the  prime  mover  stack.  Figure  3  shows  the 
prime  mover  configuration. 

Both  the  heat  exchangers  and  prime  mover  stack  have  parallel-plate  geometry.  The 
work  output  of  this  prime  mover  stack  is  modeled  by  Eq.  (18)  in  Chapter  II.  The 
thermo-viscous  losses  in  heat  exchangers  are  modeled  using  Eq.  (19).  The  losses  in  the 
cold  and  hot  ends  are  calculated  using  Eq.  (19).  Equation  (20)  and  (21)  were  used  in  the 
calculation  instead  of  Eq.  (12)  and  Eq.  (13). 

Figures  4  to  7  compare  the  MATLAB  program  results  with  the  experimental  results. 

The  experimental  data  presented  in  these  figures  are  composites  of  individual  data  sets 
recorded  at  different  times.  Figures  4  and  5  show  how  1/Q  changes  with  the  AT  for  the 
prime  mover  filled  with  helium  gas  at  a  mean  pressure  of  376  kPa  and  500  kPa.  The 

computer  simulation  tends  to  underpredict  the  experimental  onset  values  by  about  20-30 
degrees.  Figures  6  and  7  show  how  resonance  frequency  changes  with  AT.  The 
difference  between  simulation  and  experiment  is  less  than  6  Hz,  or  approximately  1%. 

One  known  limitation  of  the  standing  wave  analysis  is  that  it  ignores  boundary  layer  effect 
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Figure  3  -  Lin's  prime  mover  configuration 
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0.03 


Lin's  prime  mover  (pressure=376  kPa) 


Figure  4  -  Graph  of  1/Q  vs.  Ar  for  Lin's  prime  mover  filled 
with  helium  gas  at  a  mean  pressure  of  376  kPa. 
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Un‘s  prime  mover  (pressure=376  kPa) 


Figure  6  -  Graph  of  resonance  frequency  vs.  AT  for  Lin's  prime  mover 
filled  with  helium  gas  at  a  mean  pressure  of  376  kpa. 
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when  calculating  resonance  frequency.  However,  the  results  are  acceptable  for  the 
intended  purpose. 

B.  DELTAE  PROGRAM  SIMULATION 

The  results  of  the  DeltaE  simulation  of  Lin's  prime  mover  are  shown  in  Figs.  8-10. 
The  prime  mover  is  filled  with  helium,  at  a  mean  pressure  of  376  kpa.  Figure  8  shows  the 
heat  input  to  hot  heat  exchanger  as  a  function  of  the  work  output  of  the  prime  mover. 
Notice  that  there  is  a  linear  relationship.  Extrapolation  of  a  linear  least  squares  fit  to  the 

points  shows  that  the  work  output  goes  to  zero  at  a  heat  input  of  approximately  44  W. 
Figure  9  shows  the  square  of  the  acoustic  pressure  amplitude  as  a  function  of  .  There 
is  also  a  linear  dependence  between  the  two  variables.  Figure  10  shows  AT  as  a  function 
of  .  AT  has  a  minimum  value  of  approximately  321  K,  at  a  heat  input  of 
approximately  70  W.  This  minimum  is  taken  to  be  the  ATonset-  This  value  is  indicated  in 
Fig.  4  by  the  "  *  "  symbol.  The  DeltaE  result  is  in  good  agreement  with  both  the 

measured  result  and  the  result  of  the  MATLAB  program. 

A  similar  procedure  (see  Figs.  1 1- 13)  yields  a  ATonset  of  approximately  341  K  for  a 
mean  gas  pressure  of  500  kPa.  This  result  is  indicated  on  Fig.  5  by  the "  *  "  symbol.  This 

value  is  also  in  good  agreement  with  the  measured  predicted  results.  Detailed  input  files 
for  DeltaE  are  given  in  Appendix  A  and  B. 
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QH 

Figure  9-  Graph  of  the  square  of  the  acoustic  pressure  amplitude 
vs.  heat  input  for  Lin's  prime  mover  filled  with  helium 
gas.  The  results  are  from  DeltaE.  The  mean  gas  pressure 
is  376  kPa. 
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Figure  1 1  -  Graph  of  heat  input  vs.  work  output  for  Lin's  prime  mover 
filled  with  helium  gas.  The  results  are  from  DeltaE.  The 
mean  gas  pressure  is  500  kPa. 
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Figure  12  -  Graph  of  the  square  of  the  acoustic  pressure  amplitude 
vs.  heat  input  for  Lin's  prime  mover  filled  with  helium  gas 
The  results  are  from  DeltaE.  The  mean  gas  pressure  is 
500  kPa 
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IV.  TOROIDAL  PRIME  MOVER 


Having  demonstrated  the  ability  to  predict  the  onset  conditions  for  Lin's  prime 
mover  with  the  MATLAB  program,  the  next  step  is  to  analyze  a  toroidal  prime  mover. 
One  possible  configuration  is  shown  in  Fig.  14.  Heat  is  supplied  to  the  hot  heat  exchanger 
and  heat  is  removed  from  the  cold  heat  exchanger,  forcing  a  temperature  gradient  across 
the  stack  as  usual.  A  rigid-ended  prime  mover  would  simply  have  rigidly-capped  hot  and 
cold  (isothermal)  ducts  attached  to  the  heat  exchanger.  However,  in  the  toroidal 
geometry  the  hot  and  cold  ducts  are  connected  together,  resulting  in  a  temperature 
discontinuity  at  the  junction.  Of  course,  this  discontinuity  would  be  smoothed  out  by 
thermal  conduction,  so  that  there  would  be  a  gradual,  continuous  temperature  change 
along  the  duct.  This  gradient  will  have  some  thermoacoustic  effect  that  must  be  taken  into 
account.  More  importantly,  if  the  duct  has  a  non-negligible  thermal  conductivity,  a 
significant  amount  of  heat  will  be  conducted  through  it,  perhaps  requiring  an  excessive 
heat  input  to  reach  onset. 

To  avoid  the  detrimental  conduction  through  the  duct,  a  toroidal  prime  mover  might 
have  two  insulator  sections,  one  attached  to  each  of  the  hot  and  cold  heat  exchangers. 
The  use  of  insulators  will  allow  the  duct  to  be  held  at  ambient  temperature  while  at  the 
same  time  having  little  conduction  through  it.  The  disadvantage  of  the  insulators  is  that 
the  temperature  gradient  across  them  can  be  high.  However,  since  the  surface  area  is 
much  smaller  than  that  of  the  stack,  any  thermoacoustic  effect  will  be  small. 

Due  to  unfortunate  time  constraints,  attemps  to  build  a  toroidal  prime  mover  were 
limited  to  ofif-the-shelf  components  and  those  already  on  hand.  This  combination  resulted 
in  a  compromised  design.  These  compromises  were  understood  from  the  outset,  yet  it 
was  decided  to  complete  the  construction  as  a  learning  step  toward  a  fimctioning  toroidal 
prime  mover. 

This  chapter  contains  a  description  of  the  toroidal  prime  mover,  analysis  of  the 
prime  mover  based  on  the  MATLAB  program,  and  a  description  of  measurements 
performed,  along  with  the  results. 
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Figure  14  -  Toroidal  Prime  Mover  Configuration 
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A.  TOROroAL  PRIME  MOVER  DESIGN 


Figure  14  shows  an  illustration  of  the  toroidal  prime  mover.  This  section  contains 
specific  details.  The  hot  and  cold  heat  exchangers  and  the  stack  and  stack  holder  are 

those  used  by  Castro  [Ref  2],  More  complete  details  can  be  found  there. 

1.  Ambient  Duct  Assembly 

The  ambient  duct  assembly  is  made  from  four  sections  of  bent,  galvanized  steel 
conduit.  Each  section  is  4.4  cm  O.D.,  4  cm  I.D.,  and  48  cm  long.  The  tubes  are  joined 

together  with  2.35  cm  long  brass  collars  and  soldered. 

2 .  Insulator 

The  stainless  steel  insulators  are  sections  of  tube,  4  cm  long,  3.85  cm  I.D.  and  4.0 
cm  O.D.  They  are  located  between  the  ambient  duct  assembly  and  hot  and  cold  heat 
exchangers  to  reduce  heat  conduction  to  the  duct.  The  stainless  steel  tube  is  soldered  to 
two  stainless  steel  flanges.  One  flange  is  bolted  to  the  ambient  duct,  the  other  to  the  heat 

exchanger  and  stack  holder. 

3.  Hot  Heat  Exchanger  Assembly 

The  assembly  consists  of  hot  heat  exchanger,  hot  heat  exchanger  holder  and  hot  heat 
exchanger  heater  collar.  The  copper  heat  exchanger  is  3.924  cm  in  diameter  and  0.257  cm 
thick.  Each  plate  within  the  heat  exchanger  is  0.0292  cm  thick.  The  plate  gap  is  0.0584 
cm.  The  heat  exchangers  are  soldered  into  the  heat  exchanger  holder.  The  copper  heat 
exchanger  holder  is  7.366  cm  O.D.  and  0.953  cm  thick.  There  is  an  0-ring  groove  in  the 
side  facing  the  stack. 

The  hot  heat  exchanger  heater  collar  is  made  of  0.953  cm  thick  copper  and  bolted 
around  the  hot  heat  exchanger  holder.  The  collar  is  roughly  square  in  cross  section,  8.2 
cm  long  on  edge.  It  has  a  7.37  cm  I  D.  hole  cut  in  it  to  fit  over  the  heat  exchanger  holder. 
The  collar  is  split  and  bolted  together  around  the  holder.  There  are  two  0.625  cm 
diameter  holes  drilled  in  each  piece  to  hold  two  heater  cartridges,  used  to  maintain  a 

constant  hot  heat  exchanger  temperature. 

4.  stack  Assembly 

The  prime  mover  stack  assembly  consists  of  the  stack  and  the  stack  holder.  The 
stack  is  made  of  0.01  cm  thick  polyester  film.  The  film  is  rolled  spirally  around  a  phenolic 
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core  0.635  cm  in  diameter  and  2.529  cm  long.  The  polyester  film  is  separated  by  0.851 

mm  diameter  fishing  line.  The  stack  is  2.553  cm  long,  3,772  cm  O.D.  and  has  a  total 
surface  area  of  583  cm^ . 

The  304  stainless  steel  stack  holder  is  2.576  cm  long  and  3.772  cm  I.D.  (The  stack 
holder  is  0.22  mm  longer  than  the  stack  creating  a  0. 1 1  mm  gap  on  each  side  of  the  stack 
inside  the  holder).  The  stack  holder  has  eight  8-32  machine  screw  holes  in  each  flange  to 
attach  the  insulator  and  hot  heat  exchanger  on  one  side  and  the  cold  heat  exchanger  and 
insulator  on  the  other. 

5. Cold  Heat  Exchanger  Assembly 

The  assembly  consists  of  cold  heat  exchanger,  cold  heat  exchanger  holder  and  cold 
finger.  In  order  to  impose  a  temperature  gradient  across  the  stack,  a  65.5  cm  long,  0.7  cm 
thick,  10  cm  wide  copper  plate  (or  cold  finger)  is  bolted  around  the  cold  heat  exchanger 
holder  and  immersed  in  liquid  nitrogen.  The  cold  heat  exchanger  assembly  is  similar  to  the 
hot  heat  exchanger  assembly,  except  that  the  cold  heat  exchanger  holder  has  no  0-ring 
groove.  The  cold  heat  exchanger  is  sealed  by  using  a  0.03  cm  thick  teflon  gasket  in  the 
stack  side  of  the  cold  heat  exchanger  holder  and  a  0.17  cm  diameter  lead  wire  diameter 
gasket  on  the  insulator  side. 

B.  MATLAB  PROGRAM  SIMULATION 

In  order  to  determine  the  best  operating  mean  pressure,  resulting  in  the  lowest 
onset,  the  toroidal  prime  mover  has  modeled  with  the  MATLAB  program.  As  was 
mentioned  earlier,  designing  a  functioning  toroidal  prime  mover  is  considerably  more 

difficult  than  designing  a  conventional  prime  mover.  For  example,  in  a  conventional  prime 
mover,  placement  of  a  rigid  boundary  in  the  engine  insures  that  the  acoustic  velocity  is 
zero  at  that  point.  These  boundary  conditions  serve  as  a  convenient  starting  point  for 
computational  analysis  of  the  performance  of  an  engine.  Unlike  the  conventional  prime 
mover,  the  toroidal  prime  mover  has  no  rigid  boundary.  Therefore,  the  location  of  the 
stack  relative  to  the  nearest  pressure  antinode  is  an  adjustable  parameter. 

Figure  15  shows  a  3-D  plot  of  A  render  versus  mean  gas  pressure  of  neon  and  Xs  , 
where  Xs  is  the  distance  from  the  location  of  the  assumed  pressure  antinode  to  the  hot  end 
of  the  stack.  The  shape  of  the  surface  indicates  the  general  dependence  of  txT onset  on 
those  two  variables.  It  is  seen  that  the  value  of  jCj  yielding  the  lowest  hT onset  is  21  cm. 
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onset  temperature  difterence(K) 


207  cm  toroidal  prime  mover(pms  length=2.529cm  rad=1 .905  cm) 


Figure  15  -  Graph  of  Monset  vs.  mean  pressure  and  stack  location. 
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Figure  16  shows  ATonset  versus  mean  gas  pressure  forxj  =21  cm.  It  is  seen  that  the 
lowest  value  of  ATo„set  is  approximately  190  K  at  a  pressure  of  100  kPa.  Acknowledging 
the  tendency  of  the  program  to  under  predict  ATomet,  a  reasonable  value  for  ATonset  is 
approximately  190-210  K.  The  largest  possible  AT  attainable  with  liquid  nitrogen  and  a 

hot  end  at  nominal  room  temperature  is  approximately  220  K.  Therefore,  it  is  seen  that 
the  expected  onset  conditions  are  right  at  the  limits  of  what  is  easily  achieved. 

C.  EXPERIMENTAL  PROCEDURE 

The  toroidal  prime  mover  was  pumped  down  using  a  vacumn  pump  to 
approximately  10  kPa  and  filled  back  with  neon  to  approximately  100  kPa.  This 
procedure  is  performed  three  times  to  ensure  the  gas  inside  the  prime  mover  is 
substantially  pure  neon.  Next,  the  dewar  was  filled  with  liquid  nitrogen,  cooling  the  cold 
heat  exchanger.  As  the  cold  heat  exchanger  cooled  down,  heat  was  supplied  to  the  hot 
heat  exchanger  via  the  heater  cartridges  to  maintain  the  hot  end  temperature  at 
approximately  293  K.  The  cold  heat  exchanger  was  surrounded  by  thermal  insulation  to 
reduce  the  heat  leak  to  the  surroundings.  The  temperature  of  the  cold  heat  exchanger  was 
measured  with  a  type  J  thermocouple  clamped  between  the  cold  heat  exchanger  holder 
and  the  cold  finger  collar.  The  coldest  temperature  achieved  was  approximately  1 10  K. 

Temperature  control  of  the  hot  heat  exchanger  is  achieved  by  a  digital  temperature 
controller  (Omega  CN  9000  temperature  controller),  heater  cartridges,  and  a  power 
supply  amplifier  (Techron  5530  power  supply  amplifier).  The  reference  temperature  for 
the  temperature  controller  is  taken  from  a  type  J  thermocouple  placed  between  the  hot 
heat  exchanger  holder  and  the  heater  collar.  The  controller  sends  a  10  volt  DC  pulse  to 
the  mono  input  in  the  power  supply  amplifier.  The  duty  cycle  of  the  phase  is  proportional 
to  the  difference  between  the  set  point  temperature  and  the  thermocouple  reading.  Power 

to  the  heater  catridges  is  provided  by  the  output  of  the  power  supply  amplifier.  The  hot 
end  temperature  was  maintained  to  within  approximately  ±5  K  of  the  set  point 
temperature,  initially  293  K. 

The  pressure  inside  the  toroidal  prime  mover  is  measured  with  a  piezoresistive 
transducer  with  a  sensitivity  of  17.05  mv/psi.  The  transducer  is  mounted  in  the  wail  of  the 
ambient  duct,  approximately  40  cm  from  the  hot  end  of  the  stack.  The  transducer 
measures  the  mean  pressure  as  a  DC  signal  and  the  acoustic  pressure  amplitude  as  an  AC 
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207  cm  toroidal  prime  mover  prime  mover 


Figure  16  -  Graph  of  onset  vs.  mean  pressure  of  neon  gas  for  Xj  =  21  cm 
as  predicted  with  the  MATLAB  program. 


signal,  simultaneously.  The  transducer  output  is  sent  to  a  DC  coupled  preamplifier  with  a 
gain  of  10. 

D.  EXPERIMENTAL  RESULTS 

The  toroidal  prime  mover  failed  to  reach  onset  at  a  mean  pressure  of  100  kPa  and 
a  AT  of  approximately  180  K  (cold  end  temperature  =  1 10  K,  hot  end  temperature  =293 
K).  At  this  point  the  hot  end  temperature  was  increased  to  approximately  330  K,  and  the 

mean  pressure  varied  over  a  range  from  approximately  20  kPa  to  140  kPa.  The  prime 
mover  still  failed  to  go  into  onset. 

There  are  two  possible  explanations.  1)  The  applied  ^T  was  not  large  enough,  or  2) 
The  prime  mover  would  never  go  into  onset  for  reasons  discussed  earlier  dealing  with  how 
the  stack  might  influence  the  orientation  of  the  standing  wave. 

The  null  result  suggests  two  possible  paths  for  further  study.  1)  Design  a  prime 
mover  with  a  lower  ^Tonset  ■  This  path  is  the  subject  of  Chapter  V.  2)  Modify  the  existing 
toroidal  prime  mover  to  make  it  go  into  onset,  to  test  the  computer  model. 

It  was  decided  to  remove  approximately  half  the  ambient  duct  from  the  toroidal 
pnme  mover  and  rigidly  cap  the  remaining  sections.  This  operation  converted  the  toroid 
to  a  half-toroid.  The  new  prime  mover  configuration  is  similar  to  a  conventional  prime 
mover  except  for  the  insulator  sections  and  the  two  end  ducts  having  the  same 
temperature. 

The  duct  was  cut  and  capped  such  that  the  hot  end  of  the  stack  was  20  cm  from  the 

rigid  end.  The  total  length  of  the  half-toroid  was  108  cm.  Figure  17  shows  the  output  of 
the  MATLAB  program  for  the  half-toroid  configuration.  The  predicted  ATonsei  is 
approximately  1 17  K  at  a  mean  pressure  of  100  kPa. 

The  prime  mover  was  filled  with  neon  gas  to  a  pressure  of  100  kPa  and  the  cold 

end  was  cooled  down  with  liquid  nitrogen,  as  explained  previously.  Onset  was  observed 
at  a  A  r  of  approximately  1 10  K.  This  measured  result  agrees  with  the  predicted  value  to 
within  7  K.  The  agreement  gives  some  indication  of  how  much  the  prediction  for  the 
whole-toroid  might  be  in  error. 
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temperature  difference  (K) 


half-toroid  prime  mover 


Pressure  (kPa) 

Figure  17  -  Graph  of  Aro„,e,  vs.  mean  pressure  of  neon  gas  for  the 
half-toroid  prime  mover  as  predicted  with  the  MATLAB 

program.  The  "o"  represents  the  measured  onset  point. 
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V.  IMPROVED  DESIGN 


A.  OPTIMIZING  PARAMETERS 

The  null  result  of  the  attempt  to  reach  onset  with  the  toroidal  prime  mover  discussed 

in  Chapter  IV,  leaves  the  basic  question  "Will  a  toroidal  prime  mover  work?"  unanswered. 
One  obvious  step  is  to  redesign  the  prime  mover  to  lower  its  predicted  ATomet  to  a  value 
easily  attainable.  This  redesign  is  the  topic  of  this  chapter. 

In  considering  the  design,  it  was  decided  to  assume  the  same  stack  and  heat 
exchanger  plate  spacing  as  described  in  Chapter  IV.  The  heat  exchanger  length  is  also  the 

same.  Three  parameters  were  investigated:  1)  the  length  of  the  stack,  2)  the  radius  of  the 

duct,  and  3)  the  circumference  of  the  toroid.  It  should  be  pointed  out  that  ATonset  for 
some  of  the  parameters  shown  exceeds  216  K  (the  value  of  AT  possible  between  nominal 

room  temperature,  293  K,  and  liquid  nitrogen  temperature,  77  K).  Therefore,  for  these 
plots  we  have  set  the  hot  heat  exchanger  temperature  equal  to  373  K.  The  temperature  of 
the  ambient  duct  is  is  293  K.  Thus,  there  is  a  100  K  temperature  difference  across  the  hot 

end  insulator.  This  temperature  difference  is  in  the  wrong  direction,  thermoacoustically. 
This  will  tend  to  make  the  predicted  ATonset  higher  than  might  be  expected.  Also,  the 
thermophysical  parameters  of  the  gas  will  be  somewhat  different.  However,  these 

considerations  do  not  significantly  alter  the  conclusions  draw  from  the  results. 

Figure  1 8  shows  the  effect  of  changing  the  length  of  the  stack.  The  graph  shows 
ATonset  versus  mean  pressure  of  neon  for  Xs  =  21  cm.  The  symbols  are  for  three  different 
stack  lengths:  3.175  cm,  2.529  cm,  and  1.905  cm.  The  circumference  is  207  cm.  It  is 

seen  that  the  shorter  stack  gives  the  lowest  value  of  ATonset-  There  are  two  reasons  for 
this  result.  First,  for  a  fixed  AT,  a  shorter  stack  results  in  a  larger  temperature  gradient. 

Also,  the  shorter  stack  has  less  viscous  loss. 

Figure  19  shows  the  effect  of  changing  the  radius  of  the  prime  mover.  Again,  the 

graph  shows  ATonset  versus  gas  pressure  at  x,  =  21cm  and  a  1.905  cm  long  stack.  Three 
different  radii  are  shown:  2.54  cm,  1.905  cm,  and  1.27  cm.  The  lowest  ATonset  is  occurs 

for  the  largest  radius  and  is  approximately  133  K.  The  reason  for  the  decrease  with 

increasing  radius  is  that  1/Q  is  proportional  to  the  ratio  of  the  volume  of  gas  contributing 
to  losses  to  the  volume  of  gas  storing  energy.  This  ratio  is  roughly  proportional  to  5/7?. 
Finally,  the  effect  of  changing  the  circumference  is  investigated.  As  discussed  in 

Chapter  I,  any  nonuniformities  in  the  cross  section  of  the  toroid  may  affect  the  spatial 
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orientation  of  the  standing  wave.  Therefore,  it  is  important  to  have  as  nearly  a  uniform 
cross  section  as  possible.  One  way  to  make  a  uniform  cross  section  is  to  machine  a  toroid 
out  of  a  single  piece  of  stock.  A  207  cm  circumference  is  perhaps  too  large  to  do  this 
easily.  Therefore,  the  effect  of  decreasing  the  circumference  was  investigated.  In 
decreasing  the  circumference  from  207  cm  to  1 50  cm,  the  resonance  frequency  increases. 

Therefore,  the  position  of  the  stack  relative  to  the  standing  wave  changes.  In  order  to 

keep  the  stack  at  the  right  place,  the  value  of  Xj  must  be  decreased  by  a  factor  150/207  to 
keep  kxs  the  same  in  the  two  cases.  Therefore,  Xs  for  the  150  cm  circumference  resonator 

is  15  cm.  The  results  of  this  simple  calculation  are  verified  in  Fig.  20  which  is  a  3-D  plot 

for  the  150  cm  prime  mover.  It  is  seen  that  the  minimum  value  of  ATonset  (152  K)  occurs 
at  Xs  approximately  equal  to  15  cm.  Figure  21  shows  ATonset  versus  gas  pressure  for  a 

1.905  cm  long  stack  and  2.54  cm  radius  toroid.  Two  circumferences  are  shown;  207  cm 
and  150  cm.  It  is  seen  that  decreasing  the  circumference  increase  ATonset  to  approximately 
152  K.  It  also  decreases  the  optimum  mean  pressure  to  approximately  70  kPa. 

Decreasing  the  circumference  more  increases  ATonset  further.  The  150  cm  circumference 
is  a  reasonable  compromise  between  the  need  for  a  low  ATonset  and  having  a  prime  mover 

size  small  enough  for  machining  that  is  easily  from  a  single  piece  of  stock. 

B.  SUMMARY 

The  main  finding  of  this  thesis  is  a  preliminary  design  of  a  toroidal  prime  mover. 

The  results  of  a  first-step  design  indicates  that  a  useful  geometry  has  a  1.9  cm  long  stack 
and  2.5  cm  radius.  ATonset  for  this  configuration  is  predicted  to  be  approximately  150  K  at 
70  kPa  mean  pressure  of  neon. 
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The  effect  of  changing  the  length  of  stack 
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Figure  1 8  -  Graph  of  ATo„set  vs.  mean  gas  pressure  for  the  toroidal 
prime  mover  for  three  stack  lengths. 


The  effect  of  changing  the  radius  of  stack 
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Figure  19  -  Graph  of  onset  vs.  mean  gas  pressure  for  the  toroidal 
prime  mover  for  three  different  duct  radii. 
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Figure  21  -  Graph  of  ATonset  vs.  mean  gas  pressure  for  the  toroidal 
prime  movers  for  two  different  circumferences. 


APPENDIX  A.  INPUT  FILE  FOR  LIN'S  PRIME  MOVER 


TITLE  Lin's  Thermoacoustic  Prime  Mover  Input  File 
BEGIN  Initial  0 

3.76E+05  a  Mean  P  Pa 
591.  bFreq.  Hz  G 
400  c  T-beg  K 
3.76E+04  d  lpl(a)0  Pa  G 
0.000  e  Ph(p)0  deg 
0.000  fiUI@0  m^3/s 
0.000  g  Ph(U)0  deg 

helium  Gas  type 
copper  Solid  type 

ENDCAP  Hot  End  1 

1 . 1 46E-03  a  Area  m^2 

sameas  0  Gas  type 
copper  Solid  type 

ISODUCT  Hot  Duct  2 
sameas  la  a  Area  m'''2 
0,12  b  Perim  m 
0.055  c  Length  m 

sameas  0  Gas  type 
copper  Solid  type 

HXFRST  HotHX  3 

sameas  la  a  Area  m^2 
0.697  b  GasA/A 
7.620E-03  c  Length  m 
5.200E-04  d  yO  m 
250.0  e  Heatin  W 
500.0  fEst-T  K  (t) 
sameas  0  Gas  type 
nickel  Solid  type 

STKSLAB  Prime  Mover  Stack  4 
sameas  la  a  Area  m'^2 
0.733  b  GasA/A 
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3.500E-02  c  Length  m 
3.850E-04  d  yO  m 
1 .400E-04  e  Lplate  m 
helium  Gas  type 
stainless  Solid  type 

HXLAST  Ambient  HX 
sameas  la  a  Area  m'^2 
0.697  bGasA/A 
2.040E-02  c  Length  m 
5.2E-04  d  yO  m 
0.000  e  Heatin  W  (t) 
293.  fEst-T  K  (t) 
helium  Gas  type 
copper  Solid  type 

ISODUCT  Cold  Duct 


sameas 

5a  a  Area  m^2 

0.12 

b  Perim  m 

0.8795 

c  Length  m 

helium 

Gas  type 

ideal 

Solid  type 

ENDCAP  Cold  End 

sameas 

5a  a  Area  m'^2 

helium 

Gas  type 

copper 

Solid  type 

HARDEND  8 

0.000 

a  R(l/Z) 

0.000 

bI(l/Z) 

helium 

Gas  type 

copper 

Solid  type 

APPENDIX  B.  INPUT  FILE  FOR  LIN'S  PRIME  MOVER 


TITLE  Lin's  Thermoacoustic  Prime  Mover  Input  File 
BEGIN  Initial  0 

5E+05aMeanP  Pa 
591.  bFreq.  Hz  G 
400  c  T-beg  K 
5E+04  d  lpl(a)0  Pa  G 
0.000  e  Ph(p)0  deg 
0.000  flUI@0  m'^3/s 
0.000  g  Ph(U)0  deg 
helium  Gas  type 
copper  Solid  type 

ENDCAP  Hot  End  1 

1 . 1 46E-03  a  Area  m'^2 

sameas  0  Gas  type 
copper  Solid  type 

ISODUCT  Hot  Duct  2 

sameas  la  a  Area  m^2 
0.12  bPerim  m 
0.055  c  Length  m 

sameas  0  Gas  type 
copper  Solid  type 

HXFRST  HotHX  3 

sameas  la  a  Area  m'^2 
0.697  b  GasA/A 
7.620E-03  c  Length  m 
5.200E-04dyO  m 
250.0  eHeatIn  W 
500.0  fEst-T  K  (t) 
sameas  0  Gas  type 
nickel  Solid  type 

STKSLAB  Prime  Mover  Stack  4 
sameas  la  a  Area  m^2 
0.733  b  GasA/A 
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3.5(X)E-02  c  Length  m 
3.850E-04dy0  m 
1 .400E-04  e  Lplate  m 
helium  Gas  type 
stainless  Solid  type 

HXLAST  Ambient  HX 
sameas  la  a  Area  m'^2 
0.697  bGasA/A 
2.040E-02  c  Length  m 
5.2E-04dy0  m 


0.000 

e  Heatin  W  (t) 

293. 

fEst-T  K  (t) 

helium 

Gas  type 

copper 

Solid  type 

ISODUCT  Cold  Duct 

sameas 

5a  a  Area  m'^2 

0.12 

b  Perim  m 

0.8795 

c  Length  m 

helium 

Gas  type 

ideal 

Solid  type 

ENDCAP  Cold  End 

sameas 

5a  a  Area  m''2 

helium 

Gas  type 

copper 

Solid  type 

HARDEND  8  8 

0.000 

a  R(l/Z) 

0.000 

bI(l/Z) 

helium 

Gas  type 

copper 

Solid  type 
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